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Abstract 

Microscopic theories beyond mean-field are developed to include pairing, 
in-medium nucleon-nucleon collisions as well as effects of initial fluctuations 
of one-body observables on nuclear dynamics. These theories are applied 
■ to nuclear reactions. The role of pairing on nuclear break-up is discussed. 

. By including the effect of zero point motion of collective variables through a 

stochastic mean-field theory, not only average evolution of one-body observ- 
ables are properly described but also fluctuations. Diffusion coefficients in 
fusion as well as mass distributions in transfer reactions are estimated. 



1 Introduction 

Time-Dependent Energy Density Functional (TD-EDF) provides a suitable micro- 
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^ \ scopic framework to treat both nuclear structure and reactions along the nuclear 

^ ' chart [1,2] . In most of current applications, an effective interaction (of Skyrme or 

Gogny type) is introduced to provide an energy functional, denoted £{p)^ where p 
^ ' is the one-body density matrix. Then, guided by the Hamiltonian case, equations 

^ . of motion are written in terms of the one-body density evolution 
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where h[p\ = dS{p)/dp denotes the mean-field Hamiltonian. Although this approach 
O ■ is traditionally called Time-Dependent Hartree-Fock (TDHF), it should not be con- 

fused with TDHF derived from a two-body Hamiltonian. At least because of the 
parameters of the effective vertex are directly adjusted to experimental observations. 



' EDF incorporates much more correlations than a pure Hartree-Fock theory. TD- 

EDF has been improved significantly in the past decades and could now be applied 
without assuming specific symmetries in space and includes all terms of the effective 
interactions used in the static EDF [1]. 

Starting from a pure independent particle state, the one-body evolution described 
by Eq. ([T]) can be replaced by the evolution of a single Slater determinant 0. As a 



^This is due to the fact that — p is preserved along the mean- field trajectory. 



consequence it misses important physical effects which are accounted for in "state 
of the art" EDF dedicated to nuclear structure. These effects are (i) Pairing cor- 
relations usually treated by considering quasi-particles trial states instead of Slater 
determinants, (ii) Fluctuations of one-body observables and correlations induced 
by restoration of broken symmetries which are generally incorporated using config- 
uration mixing methods, leading to the so-called Multi-Reference EDF (MR-EDF) 
[3-5]. Due to the numerical effort, Time-Dependent Energy Density Functional de- 
voted to nuclear dynamics generally neglects these correlations and fails to account 
for the richness of phenomena taking place in nuclear dynamics [6]. Recently, we 
have developed two transport theories dedicated to effect (i) and (ii) respectively. 
Highlights and applications of these new approaches are given below. 



2 Time-Dependent EDF with pairing correlations 

Guided by the Hamiltonian case, different extensions of mean-field theory have been 
proposed starting from a generalization of the one-body density equation of motion: 

thf^=[h[p],p] + lTT,[v'i„C,,], (2) 

where denotes the anti-symmetric effective vertex in the correlation channel. 
Tr2(.) is the partial trace on the second particle while Cu = pu — PiP2(l — ^12) 
denotes the two-body correlation matrix defined from the two-body density pi2. 
Here, indices refer to the particle on which the operator is applied while A12 is the 
permutation operator. Eq. should be complemented by the equation of motion 
for the two-body correlation. It is written here as 

'^~d^ = [h[p]+h2[p],C,2] 

+ (1 - Pl)(l - P2)^j2PlP2 - PlP2^^?2(l - Pl)(l - P2) ^ B12 
+ (1- pi- p2)Vi2Cl2-Ci2Vi2i'^- Pl- P2)- ^ Pl2 

+ Higher orders. (3) 

B12 and P12 contain correlations associated to in-medium nucleon-nucleon collisions 
and pairing respectively. Higher-order terms (not shown here) contain in particular 
three-, four-... body effects. When correlations between more than two particles 
are neglected, the resulting theory is known as the Time-Dependent Density-Matrix 
(TDDM) (see for instance [7]). Even in that case, prohibitory numerical efforts are 
needed due to the explicit treatment of two-body matrices. Guided by the BCS 
approximation, we recently proposed to reduce the complexity by assuming that 
components of VI2 and Cu are non-zero only between pairs of time-reversed states. 
This approximation, called hereafter TDDM^, leads to important simplifications: 
(a) the number of correlation matrix elements to be calculated is significantly re- 
duced, and (b) H12 cancels out and only B12 and Pu contribute to the correlation 



evolution. The TDDM^ theory has been recently applied to the nuclear break-up of 
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Figure 1: (Color online) Left: One body density for different times of the dynamical 
evolution for an ^^O +^°®Pb calculation at 40 A.MeV. The circle represents the 
^'^^Pb target. Right: Final distribution of relative angles between emitted nucleons 
initially correlated (black) or anti-correlated (red). (Adapted from ref. [8]) 



correlated systems. Figure [T] illustrates the evolution of the one-body density in a 
^^O as a ^ospb collision partner passes by. The perturbation by the external nuclear 
passing potential induces an emission to the continuum. The one nucleon emission 
has been extensively studied in ref. [9]. With the TDDM-^ theory, one can extend 
the study of one nucleon case to the two nucleons emission for initially correlated 
nucleons. In left side of figure [H the relative angle between two nucleons emitted 
in coincidence has been extracted for initially correlated or anti-correlated systems. 
As expected intuitively [8,10], in the former case small relative angles are obtained 
while larger relative angles are seen in the latter case. Our study clearly points out 
that coincidence measurements using nuclear break-up can be used as a tool to infer 
correlation properties inside nuclei. 



3 Proper treatment of one-body observables fluc- 
tuations 

Besides the absence of pairing correlation discussed in previous section, mean-field 
transport theory given by Eq. [T] provides a good description of mean values of one- 



body observables in low energy reactions. However, it completely fails to describe 
fluctuations of one-body observables. This can directly be traced back to the use 
of a single Slater determinant as a building block. Indeed, as it is well known in 
the static case, the independent particle picture is unable to describe properly the 
zero point motion in collective space pleading for multi-reference approaches where 
several Slater determinants (or quasi-particle states) are considered simultaneously. 
During the past decades, large efforts have been devoted to develop transport theo- 
ries that are able to describe not only mean- values but also fluctuations (for a review 
see [6,11]). However, even 30 years after the first application of TD-EDF, no practi- 
cal solution simple enough to be applied to many different physical phenomena has 
been really proposed so far. Recently, we have shown that the theory introduced in 
ref. [12] and based on stochastic initial values in one-body space might provide such 
a solution. The underlying principle is the following: 

Imagine a correlated system, and some one-body observable Q with mean-value {Q) 
and fluctuation gq — {Q^) — {QY- A single Slater determinant (SD) that minimizes 
the EDF is generally able to reproduce the mean- value but strongly underestimates 
fluctuations, i.e. Gq^ <^ ctq, where ctq^ denotes the Slater determinant expecta- 
tion value. If instead, a statistical ensemble of Slater determinants is used. Each 
SD leads to a quantum expectation value cxg, where A labels the SD under inter- 
est. One can then optimize the statistical ensemble in such a way that ag ~ ag 
where (X) denotes the statistical average over quantum expectation values. The 
statistical assumption described above provides a practical solution to mimic initial 
zero point motion in collective space. Then, the system evolution is performed by 
evolving each SD independently from the others, therefore neglecting interference 
between different trajectories. This approach, that can be regarded as the first step 
towards multi-reference TD-EDF, has been shown (i) to incorporates the one-body 
dissipation and associated fluctuation mechanism in accordance with the quantal 
dissipation-fluctuation relation, (ii) to give dispersion of one-body observables that 
is identical to a previous formula derived from the Balian-Veneroni variational prin- 
ciple [13,14]. 

In practice, initial fluctuations are simulated by an ensemble of initial single-particle 
density matrices, each of them associated to a single SD [12]: 



where summation i and j is made over a complete set of single-particle states 
$i(r, to), and i implicitly contains spin and isospin quantum numbers. Compo- 
nents of density matrices, pfj are time-independent random Gaussian numbers with 
mean value pfj — dijUi and a variance of the fluctuating part 5p^, is specifled by. 




(4) 



SPij^Pj'i' = T^^jj'^ii' [^i(l - rij) + nj{l - Hi)] . 



(5) 



The great advantage of the Stochastic Mean-Field (SMF) theory is that each Slater 
determinant A evolves independently from each other following the time evolution 
of its single-particle wave-functions in its self-consistent mean-field Hamiltonian, 
denoted by h{p^), according to 

ih^^^,{v,t-\) = h{p^)^,{v,t-\). (6) 

with the boundary condition $j(r, t; A) = $i(r,to)- Time evolution of mean- values 
and dispersion then express as 

Qit) = ^{Mt-A)mAt-A))pf, ^QW), <yQ = {Q\t)-Qit)Y, (7) 

where the average is made over initial conditions. 



3.1 Application to fusion and transfer reactions 

The powerfulness and applicability of the SMF have been recently illustrated in 
fusion reactions by extending the work of ref. [15,16]. Using a macroscopic reduc- 
tion of the stochastic mean-field evolution, central collisions leading to fusion have 
been mapped to a one-dimensional macroscopic Langevin evolution on the relative 
distance R between the two nuclei given by [17]: 

= -dRUiR"-) - ^{R^)R^ + ep{t)- (8) 

U{R''') and 'y{R^) denote the nuclear+coulomb potential and dissipation associated 
to one-body friction respectively and are already present at the mean-field level 
[15,16]. C,p{t) is a Gaussian random force acting on the relative motion reflecting 
stochasticity in the initial value. This fluctuating part leads to diffusion in collective 
space which can be approximated by ^p(t)^p{t') ~ 26{t — t')Dpp{R) where Dpp{R) 
denotes the momentum diffusion coefficient. The latter term is nothing but the one 
that is missing in the original theory and is of primer importance to properly describe 
observables fluctuations. An example of reduced friction P{R) = 'y{R)/fi{R) and 
diffusion coefficients Dpp{R) estimated from the macroscopic reduction of SMF is 
given in figure [2] for the head-on ^°Ca-|-^'^Ca collision. 

The possibility to estimate transport coefficients associated to fluctuation and dissi- 
pation from a fully microscopic quantum transport theory is a major breakthrough. 
However, to be really convincing, one should in addition prove that the increase of 
fluctuations is consistent with experimental observations. To prove that SMF can be 
a predictive framework, we have recently considered transfer reactions [18]. Frag- 
ment mass distributions deduced from Heavy-Ion reactions have been extensively 
studied. It is seen that the dispersion in mass scales approximately with the average 
number of exchanged nucleons. While mean-field properly describes the latter, it 



30 



y 20 

O 

S 10 



\ 


' 1 ' 


' 1 

SMF 










\ ^c.m. ~ 
1 T~---J 


100 MeV- 
^ , 1 



6000 



I 

> 



4000 



2000 



\ ' 

\ 


1 ' 1 


1 ' 

SMF 


- 


\ ^°Ca + 


4«Ca 




\ F - 

1 , 1 7 1 


100 MeV - 

L , 1 , 



10 11 



12 



R [fm] 



9 10 
R [fm] 



11 12 



Figure 2: (Color online) Evolution of reduced friction (left) and diffusion coefficient 
(right) as a function of the relative distance for the head-on ^'^Ca+'^'^Ca collision at 
center of mass energy -Ec.m. = 100 MeV. 



miserably fails to account for the dispersion. This phenomena is rather well under- 
stood in macroscopic models but has not been yet reproduced microscopically. To 
address this issue, we have considered head-on collisions below the Coulomb bar- 
rier. In that case, nuclei approach, exchange some nucleons and then re-separate. 
Similarly to the relative distance case, a macroscopic reduction onto the projectile 
(resp. target) mass, denoted by (resp. A^) can be made leading to: 

jA), = v{AU)+eA{t\ (9) 

where f (Ap, t) denotes the drift coefficient for nucleons transfer. Again, the fluctuat- 
ing term i\{t) is linked to the diffusion in mass through ^^(t)'CA(^') = 2(5(t — t')-DAA- 
Mass dispersion can then directly be estimated using: 

4^(t) ~2 f DAA{s)ds. (10) 
Jo 

An illustration of o"^^(t) for ^°Ca+^°Ca reactions is given in Fig. 13. II and compared 
to the number of exchanged nucleons, denoted by A'ex- In all cases, both quantities 
are very close from each other and lead to much higher dispersion than the original 
mean-field. Indeed, the estimated asymptotic values in the latter case are 0.004, 
0.008 and 0.008 from low to high energy and are much less than the final number of 
exchanged nucleons that are equal to 0.43, 1.44 and 3.63 respectively. On opposite, 
the predicted asymptotic mass dispersions are equal to 0.73, 1.72 and 3.79 and is 
much closer to N^x (see also figure 13.11) . This numerical test provides a strong 
support for the validity of the stochastic mean-field approach. 
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Figure 3: (Color online) Evolution of calculated in SMF approach for ^°Ca+^°Ca 
(top) at different center of mass energies. Number of exchanged particles is super- 
imposed by the filled-circles, filled-squares, and filled-triangles from high to low 
energies. (Adapted from [?]) 



4 Summary 

Two different approaches have been presented that include correlations beyond 
mean-field. These approaches differ in the strategy and in the types of correla- 
tion that are included. The TDDM^ formalism is appropriate to account for pairing 
effects in a reasonable numerical time. An illustration on break-up reactions clearly 
points out the importance of such correlations in nuclear reactions. The SMF theory 
with initial random conditions can be seen as the first step towards the inclusion of 
collective variables zero point motion. Applications to fusion and transfer reactions 
could be of particular interest to treat fluctuations around an average mean-field 
path. 
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